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Abstract 

We argue that the use of the universal unintegrated gluon distribution and the kr (or 
high energy) factorization theorem provides the natural framework for describing observ- 
ables at small x. We introduce a coupled pair of evolution equations for the unintegrated 
gluon distribution and the sea quark distribution which incorporate both the resummed 
leading ln(l/x) BFKL contributions and the resummed leading ln(Q 2 ) GLAP contribu- 
tions. We solve these unified equations in the perturbative QCD domain using simple 
parametic forms of the nonperturbative part of the integrated distributions. With only 
two (physically motivated) input parameters we find that this kx factorization approach 
gives an excellent description of the measurements of F2(x, Q 2 ) at HERA. In this way the 
unified evolution equations allow us to determine the gluon and sea quark distributions 
and, moreover, to see the x domain where the resummed ln(l/x) effects become signif- 
icant. We use kx factorization to predict the longitudinal structure function Fl(x,Q 2 ) 
and the charm component of F^ix, Q 2 )- 



1. Introduction 

The experiments at HERA have opened up the small Bjorken x regime. One of the most 
striking features of the data is the strong rise of the structure function F 2 as x decreases 
from 10~ 2 to below 10~ 4 Q. At first sight it appeared that the rise was due to the (BFKL) 
resummation of leading ln(l/a;) contributions @, However, with an appropriate choice of 
input distributions and of the starting scale for the Q 2 evolution, the observed growth can 
also be reproduced within the conventional GLAP framework which just sums up the leading 
(and next-to-leading) InQ 2 contributions. Indeed GLAP global fits exist which give a good 
description of the small x measurements of F 2 in the x range accessible at HERA j|, see 
also |6|, [7j . Is it possible to conclude that there will be no significant BFKL-type contributions 
to F 2 in the HERA small x regime or does a physically reasonable alternative description exist 
with sizeable ln(l/x) resummation contributions? Here we address this question. Clearly the 
specification of the non-perturbative input to the QCD evolution will be crucial. 

Recall that the basic dynamical quantity at low x is the gluon distribution f(x, k^) un- 
integrated over its transverse momentum kx- It is related to the conventional gluon density 
g(x,Q 2 ) by 

rQ 2 rlk 2 

xg{x,Q 2 ) = / ^ff(x,k 2 T ). (1) 

In the leading ln(l/a;) approximation f(x, k 2 -,) satisfies the BFKL equation and exhibits an x~ w ° 
growth and a diffusion in In k\ as x — > 0, where the BFKL or so-called hard Pomeron intercept 
ujq = (3as/7r)4 In 2 for fixed as- The observable quantities are computed in terms of / via the 
kr (or high energy) factorization prescription || |]. For example, the structure functions 
are given by 

Fi = Fi 9 ® f (2) 

where (g> denotes a convolution in transverse, as well as longitudinal, momentum. Here F? 9 are 
the off-shell gluon structure functions, which at lowest order are determined by the quark box 
(and crossed-box) contributions to photon-gluon fusion, see Fig. 1. 

The BFKL gluon f(x, kj.) and kr (or high energy) factorization theorem were used [|l(| 



to predict the small x behaviour of F2 prior to the measurements at HERA. The method 
used a starting distribution f(x ,k^) at, say, x = 0.01 deduced from an integrated gluon 
g(x ,Q 2 ) which had been determined in a global parton analysis of fixed-target deep inelastic 
and related scattering data. With this input the BFKL equation was solved to determine 
f(x, kj) in the small x domain (x < Xq). A major uncertainty in this procedure to predict F2 
is the treatment of the infrared region, k\ < k 2 . A recent study [T1J using this approach finds 



that the BFKL predictions for F 2 increase too steeply with decreasing x in comparison with the 
HERA measurements, and concludes that there is no evidence for the effects of the resummation 
of ln(l/x) terms. Before we accept such a conclusion we should note the limitations of this 
form of the test of the BFKL A>r-factorization approach. 
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It is assumed [11, 12 1 that in the infrared region f(x, fcf,) has the form 

2 2 k 2 + k 2 2 

f(x,k T <k ) = ^ f(x, k ) (3) 

with k 2 = 1 GeV 2 say, and where k 2 is an adjustable parameter. It turns out that there is a 
sizeable contribution from the infrared non-perturbative region which, since it is tied to f(x, k^), 
is forced to have the BFKL growth with decreasing x. The formalism does not include the GLAP 
leading InQ 2 resummations which go beyond the leading ln(l/x) approximation. Finally it 
should be noted that the BFKL equation for f(x,k^) only resums the leading order ln(l/x) 
terms. Sub-leading ln(l/a;) effects are expected to reduce the growth of / with decreasing x. 

Clearly this simplified procedure provides only a crude test of the underlying dynamics in 
the small x (HERA) domain. The main deficiencies are the treatment of the infrared region and 
the lack of a unified approach which incorporates both the BFKL ln(l/a;) and the GLAP ln(<5 2 ) 
resummations. An important development is the demonstration by || |9|, [K| that at the leading 
twist level the BFKL kx factorization approach can be reduced to the conventional collinear 
GLAP factorization in which the anomalous dimensions and coefficient functions are extended 
to include the full resummation of leading ln(l/x) terms. This motivated an informative study 
by Thorne He has shown (in a scheme independent way) that the inclusion of the ln(l/x) 
terms within the collinear factorization approach gives a satisfactory, and even an improved, 
description of the F 2 data. He assumes non-perturbative components to the input distributions 
for the observables at scale Qq which are 'flat' at small x, and which are the only contributions 
at the scale (A LL < lGeV 2 ) that denotes the boundary between the perturbative and non- 
perturbative regions. In this way he demands that the rise of the structure functions with 
decreasing x must entirely come from perturbative effects. In summary Thorne finds, within 
the collinear factorization framework, that the 'BFKL' ln(l/x) terms are not excluded, but 
rather are favoured, by the F 2 data. 

The agreement with the data obtained by Thorne |13[] is contrary to the conclusion of ref. [[U] 
which was based on the BFKL gluon and k^ factorization. This suggests that the application 
of the &t factorization approach was too simplistic. Here we re-examine and improve the 
determination of f(x, k^) and, via kr factorization, obtain a more realistic description of F 2 . 

The first improvement is that we study a 'unified' equation for f(x, k^) which incorporates 
BFKL and GLAP evolution on an equal footing flTJ]]. To be precise we solve a coupled pair of 
integral equations for the gluon and sea quark distributions, as well as allowing for the effects 
of the valence quarks. In this way we eliminate the problems of matching at x = Xq. A second 
improvement is a more physical treatment of the non-perturbative (or infrared) contributions 
to the BFKL equation and the kx factorization integrals. We shall see that the former can be 
specified entirely by the integrated gluon distribution at the scale Q 2 = k\ which marks the 
boundary of the perturbative and non-perturbative regions, whereas the integrals also need a 
non-perturbative component of the sea. In fact we find that an excellent description of the 
HERA measurements of F 2 is possible in terms of just two physically motivated parameters 
which fully determine these infrared contributions. 
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There is considerable merit in going back from the collinear to the kx factorization approach 
for the description of small x deep inelastic scattering. Indeed, in the reduction to collinear 
form we loose some of the physical structure which is contained in the gluon ladder and kx 
factorization. We discuss these limitations of the collinear approach below. 

First we note that the high energy (or low x) behaviour of the structure functions is driven 
by the BFKL gluon ladder coupled (through kx factorization) to the photon via the quark 
box. The corresponding Feynman diagram has a calculable perturbative contribution for all Q 2 
(except for Q 2 — > for massless quarks). The 'hard' or 'QCD' pomeron contribution generated 
by this ladder is thus present for all Q 2 . In particular there is a known perturbative contribution 
to the structure functions at Q 2 = k^ coming from configurations in which the gluon transverse 
momenta within the ladder lie in the perturbative domain k T > fcg. This contrasts with the 
collinear factorization approach in which it is contrived to describe the observables in terms of a 
purely non-perturbative contribution at some scale, say Q 2 = k 2 . The perturbative component, 
which must be present at Q 2 = k 2 ,, will be evident in the kx factorization approach that we 
introduce below. 

A more subtle limitation of the Renormalisation Group (RG) and collinear factorization 
approach at small x concerns the treatment of the running of as- We solve the BFKL equation 
with running, rather than fixed, coupling where a 5 depends on the local scale k\ along the 
ladder. This way of implementing the running of as is supported by the calculation of next- 
to-leading ln(l/x) effects |L5], |16[]. The solution of the BFKL equation with running as can be 
reduced to the conventional RG form using saddle point techniques. However, the saddle point 
approximation is not applicable for arbitrarily small values of x |T7], T3, [H| [Ti| . 



A third advantage of the k^ factorization approach is that it allows us to appropriately 
constrain the transverse momenta of the emitted gluons along the BFKL ladder. We are 
therefore able to quantify the effect of imposing this constraint. (Recall that in the usual 
application of the BFKL equation the gluon transverse momenta are taken to be unlimited.) 
We will see that the kinematic constraint largely subsumes the angular ordering constraint 
which is the basis of the 'CCFM' equation |20], ^TJ. The CCFM equation also incorporates 
both BFKL and GLAP evolution. 

Another difference is that the BFKL contribution is a sum over all twists, whereas when 
it is reduced to collinear form only the leading twist is conventionally retained. Finally kx 
factorization is much simpler to implement at small x than collinear factorization. We deal 
with dynamical quantities (namely the BFKL kernel and the structure function of the off- 
shell gluon) which can be calculated perturbatively. We calculate them to first order in as- 
Essentially the agln(l/:r) terms are effectively resummed by simply integrating over the entire 
k\ phase space allowed for the gluon ladder and the kr factorization integrals. 

In summary, the natural framework with which to describe observables at small x is the 
unintegrated gluon density f(x,k^) together with the fey factorization theorem. (Here we use 
it to calculate the observable deep inelastic scattering structure functions F 2 and Fl-) That is, 
at small x the distribution f(x,k^) is the basic, universal quantity which can be taken from 
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process to process. If we were to reduce this framework to a collinear factorized form then the 
"simple", but rich, physics structure of the gluon ladder is not fully taken into account and 
even may be distorted. It could be argued that most of the effects occur at subleading order in 
ln(l/a;) and since only the leading order is completely known, little is lost. However, the effects 
have a direct physical origin. They are clearly present and are expected to be the dominant 
corrections in a more complete analysis. 

Since we shall unify the BFKL and GLAP formalisms, the resulting equation for f(x,kj) 
is valid both for small x and large x. Moreover the region where the 'BFKL' ln(l/a;) effects 
become significant will be decided by the underlying dynamics (QCD). We use the formalism 
to fit to deep inelastic data and so we are able to quantify the importance of BFKL effects. 

The outline of the paper is as follows. In Section 2 we make modifications to the BFKL 
equation for the unintegrated gluon distribution f(x, k 2 ) which allow for the GLAP leading 
In Q 2 contributions and which enable all the non-perturbative effects to be encapsulated in an 
input distribution for the integrated gluon, xg(x, k$). In Section 3 we introduce the equation 
for the quark singlet (momentum) distribution E(x, Q 2 ), again paying particular attention to 
isolate the contribution to the non-perturbative region and to ensure that the perturbative 
terms are allowed to contribute for all Q 2 . In Section 4 we numerically solve the coupled in- 
tegral equations for / and E. The kx factorization theorem is used to calculate F 2 (x, Q 2 ) as a 
function of the two parameters that specify the input gluon (which we take to be 'flat' at small 
x). Optimum fits to the available F 2 data at small x are presented, and predictions are made 
for the charm component F 2 and for the longitudinal structure function Fl- Finally, in Section 
5 we present our conclusions. 

2. Unified BFKL and GLAP equation for the gluon 

We start from the BFKL equation for the unintegrated gluon distribution 

f(x,k 2 ) = f°\x,k 2 ) 

(4) 

^_ ,, 2 ., a f'dz fdk' 2 l f(x/z,k> 2 )-f(x/z,k 2 ) fjx/z, k 2 ) \ 

where % = 30:5/71 and k = kx,k' = k' T denote the transverse momenta of the gluons, see 
Fig. 1. The equation corresponds to the leading ln(l/x) approximation. 



2.1 From the BFKL to the unified equation 

In order to make the BFKL equation for the gluon more realistic and to extend its validity 
to cover the full range of x we make the following modifications. First, to incorporate leading 
order GLAP evolution, we add on to the right-hand side of (|j) the term 0] 

"i dz 



(5) 

where we have used (|1|). The -1 allows for the contribution which is already included in the 
BFKL equation. The inclusion of the additional term (|JJ) gives contributions to the gluon 
anomalous dimension 7 95 which are subleading in asln(l/x) but leading in as. These standard 
leading order GLAP contributions have an impact on the gluon intercept. They soften the 
small x rise of the gluon distribution and also change its overall normalisation. 



The second modification to (W) is the introduction of the kinematic constraint 22, 



k < k 2 /z (6) 

on the real gluon emission term, that is, on the integral over f(x/z,k' 2 ). The origin of the 
constraint is the requirement that the virtuality of the exchanged gluon is dominated by its 
transverse momentum, \k' 2 \ ~ k?. For clarity we have restored the subscript T in this equation. 
The constraint is another physically motivated, subleading correction in asln(i/x). 

Thirdly, we notice that the integration region over k' 2 in ([|) extends down to k' 2 = where 
we expect that non-perturbative effects will affect the behaviour of f(x, k' 2 ). We are only going 
to solve equation (|j) in the perturbative region, defined by k 2 > k 2 ,, so we only have to worry 
about the infrared contribution due to the real emission term from the interval < k' 2 < k^. 
We may rewrite this infrared contribution in the form 

, 2 f k o dk' 2 fx , /2 \ f k odk' 2 £ fx , /2 \ x fx , 2 \ 

The parameter k^{= Qq) denotes the border between the perturbative and non-perturbative 
regions. Its magnitude will be taken to be around lGeV 2 . 

Finally we must of course add to the right-hand side of (f|) the term which allows the quarks 
to contribute to the evolution of the gluon, that is 



2. J x dzP ^ E {- Z > k ) ^ 
where S is the singlet quark momentum distribution. To be explicit 

S(x, k 2 ) = ^2 x(q + q) + x(c + c) 

q=u,d,s 

(9) 

= V(x,k 2 ) + S uds (x,k 2 ) + S c (x,k 2 ) 

where V, S u d s and S c denote the valence, the light sea quark and the charm quark contributions 
respectively. We discuss the evolution equation for S(x, k 2 ) in the next section. 
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Gathering together all the above modifications, equation (H) for the gluon becomes 



f(x,k 2 ) = f^(x,k 2 ) 




(10) 

where now the driving term has the form 

/ (0) (x, k 2 ) = / (0) (x, k 2 ) + f dzP gg {zf-g (^ 2 ) . (11) 

In ( |i0|) we include the constraint k' 2 > k\ on the virtual, as well as the real, contributions 
in order to avoid spurious singularities at k 2 = k 2 ,. In the perturbative region, k 2 > k$, we 
may safely neglect the genuinely non-perturbative contribution f^(x, k 2 ) which is expected to 
decrease strongly with increasing k 2 . It is important to note that we have avoided the necessity 
to parametrize f(x,k 2 ) in the non-perturbative region. Equation (|H]) only involves f(x,k 2 ) 
in the perturbative domain, k 2 > /c 2 . The input fllTD is provided by the conventional gluon 
at the scale k 2 ,. That is the input to our 'unified BFKL + GLAP' equation is determined by 
the same distribution as in conventional GLAP evolution. The modifications to (|j) allow us 
to overcome the serious limitations discussed in the introduction. Surprisingly, we find that we 
can achieve an excellent description of all the deep inelastic data using the most economical 
parametrization of the input gluon 

xg(x,k 2 Q ) = N{l-xf. 

In particular the observed growth in F 2 (x,Q 2 ) with decreasing x is generated entirely by per- 
turbative (ln(l/x) and InQ 2 ) dynamics. 

It is easy to see how eq. ( |T0"D reduces to the conventional GLAP evolution equation for the 
gluon in the leading InQ 2 (or rather Ink 2 ) approximation. The leading lnfc 2 terms arise from 
the strongly ordered configuration, k 2 , k' 2 k 2 , for the real emission contributions and to 
the neglect of the virtual contributions. Then (|10[) becomes 



f(x,k 2 ) = 




(12) 
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where we have taken into account (|TT|) and the remarks concerning the omission of f(°>. Upon 
using ([J]) we see that (|T2| ) becomes 



d(xg(x 1 k 2 )) a s {k 2 



f 1 dz\p gg (z)-g(-,k 2 ) + P gq (z)z(-,k' 
J x L Z V Z / \ £ 



(9 In A; 2 2tt 
which is simply the GLAP evolution equation for the gluon. 



(13) 



2.2 Anomalous dimension of the gluon. 

We will solve ([10]) for the unintegrated gluon. However first we anticipate the general 
behaviour of the anomalous dimension of the gluon which will come from this equation. To do 
this we rewrite the equation in terms of the moment function 

J(u,k 2 ) = [ dxx"- l f{x,k 2 ). (14) 

We have 



o 



uj Jkl k u 





Q(k 2 - k' 2 ) + (k 2 /k' 2 Y &(k' 2 - k 2 ) 


-f(uj,k 2 ) 


\k' 2 -k 2 \ 



+ f(^k 2 



(Ak' A + k A Y 

+ a s {k 2 )P{u) Q^l^.k' 2 ) (15) 

where we have neglected, for simplicity, the contribution coming from the quarks and where 
P(oo) is the moment function of (zP gg (z)/6 — 1). The term in square brackets is due to the 
kinematic constraint. Without this constraint we would have 1 instead of (/c 2 //c' 2 ) w , and the 
two functions would simply sum to unity. 

For large k 2 the moment function behaves as 



f ^ k) ^{¥ ) (16) 

where, for illustration, we take fixed as- The quantity 7 55 is the anomalous dimension of 
the gluon. If we insert flUf ) into ( |i~5"D then we find, after some algebra, the following implicit 
equation for 7 99 

l--^(7 33 ^) - —P{u) = (17) 
where K, the double moment of the kernel in flT5|), is given by 

K{ u) = r dp r [^9(i -^+^-1)1-1 i \ (18) 

Jo p \ \p-l\ [4p 2 + l]iJ 
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If uj = then the expression in square brackets reduce to p 1 and we have the familiar BKFL 
result 

K( 1 ,lu = 0) = [2*(1) - #(7) - *(1 - 7)] (19) 
where \I/ is the logarithmic derivative of the Euler gamma function. 

It is clear that 7 S9 , which satisfies ([H]), is of the formQ 

-f gg (uj,as) = 7 BFKL (^f ) + a s P(uj) + higher order terms, (20) 
where 7 BFKL satisfies the usual equation 



1 _ ^W 7 BFKL^ = \ = Q (21) 
10 



The higher order terms include contributions which are sub leading in as/uj as well as in as- 
The anomalous dimension 7 9g has a branch point singularity in the u plane, whose position 



uj = uj (as) controls the small x behaviour of the gluon distribution. The inverse of ( |i4|) gives 

f{x,k 2 ) ~ x ~M*s)_ (22) 
The value of is obtained from the requirement that 



_ ^( 7 ,^ ) - ^P(u )\ = 0, (23) 



together with the equation 



1_^^( 7 ^ ) _ ^P( Wo ) = 0, (24) 

LUq 7 



see (|17"1) . Recall that in the leading ln(l/x) (or leading approximation (|23| ) reduces to 

dK(j 7 0)/d^ = 0, (25) 
which is satisfied when 7 = |. Thus from fl2l|) we obtain the well-known BFKL result that 

uj = a s K(i = ~,u; = 0) = a 5 41n2. (26) 



The relevant domain for solving the pair of equations (|i7|) and is < 7 < 1 and lg> > 0. In 
this region the Mellin transform of the non-singular part of the gluon splitting function satisfies 

P(u) < 0, (27) 

and moreover 

K(<y,u) < K(rf,Q). (28) 



1 Ii we were to replace ^ BFKL simply by the term which is leading order in as, that is ag/w, then the sum 
of the first two terms of (^0|) gives the conventional GLAP anomalous dimension. 
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Thus both the additional non-singular part of the GLAP splitting function (P gg — Q/z) and 
the kinematic constraint (which takes K(^y, 0) to K^jjU)) tend to reduce the magnitude of luq 
from the BFKL value shown in fl2"S|). These corrections are of course subleading in ln(l/x). 
Our numerical analysis with running as reflects this softening of the x~ u ° singular behaviour. 



2.3 The CCFM equation 

A more general treatment of the gluon ladder, which follows from the BKFL formalism 
is provided by the Catani-Ciafaloni-Fiorani-Marchesini (CCFM) equation based on angular 
ordering of the gluon emission along the chain |]20| , plj [24j- The equation embodies both the 



BKFL equation at small x and the conventional GLAP evolution at large x. The unintegrated 
gluon distribution / now acquires a dependence on an additional scale (which we may take to 
be Q 2 ) that specifies the maximal angle of gluon emission. The CCFM equation has the form 

f(x,k 2 ,Q 2 ) = f^(x,k 2 ,Q 2 ) 



+ a s I"- me(Q-qz)A R (z,k 2 ,q* 
Jx z J no 2 



cPq 
irq 



k 2 



-f(- 



[q + k)\q 2 



(29) 



where the theta function Q(Q — qz) reflects the angular ordering constraint on the emitted 
gluon. The so-called 'non-Sudakov' form-factor Ar is given by 



z,k 2 ,q 2 ) 



cxp 



--sj^J^eW-z^eie-q-) 



(30) 



Eq. ( f29|) contains only the singular 1/z term of the g — > gg splitting function (which is screened 
by the virtual corrections contained in Ar). Its generalisation to include the remaining parts of 
this vertex (as well as the quark contributions) is possible. Eq. ( p9|) has been solved numerically 
in the small x domain and the solution for f(x,k 2 ,Q 2 ) was presented in [21]. The CCFM 
equation, which is a generalisation of the BFKL equation, generates a steep a; - " type of 
behaviour^ but uq now acquires significant subleading ln(l/x) corrections which come from 
the angular ordering constraint P5 |. The constraint also introduces subleading terms in the 
anomalous dimension 



^99 



BFKL 

1 99 



— + a s hi — 

u / V uj 



+ 



so the angular ordering constraint which gives rise to the CCFM equation and the kinematic 
constraint (0) lead to similar effects - both give subleading ln(l/x) corrections to the "QCD 
Pomeron intercept" u;o and to the gluon anomalous dimension 7 99 . We found that the kinematic 
constraint overrides the angular ordering constraint except possibly in the large x domain when 
Q 2 < k 2 |p2| , see also f23|. Thus in our formulation we neglect the angular ordering constraint 
altogether. 



3. The equation for the quark distribution 

2 The effect on F2 is considered in p5| . 
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At small x the gluon drives the sea quark (momentum) distribution S via the g — > qq 
transition, see Fig. 1. We evaluate the effect using the hr factorization theorem. To be precise 
we use the hr factorization prescription to calculate observables (such as F 2 ) directly from 
the unintegrated gluon distribution f(x,k T ). For F 2 we interpret the result in terms of the 
sea quark distributions, implicity assuming the DIS scheme. The total sea is the sum of the 
individual quark contributions 

S(x,Q 2 ) = £S,(x,Q 2 ). 
i 

At small x the factorization theorem gives 




where Sb ox describes the quark box (and crossed box) contribution shown in Fig. 1. S^ ox 
implicitly includes an integration over the transverse momentum, k, of the exchanged quark. 
Indeed, evaluating the box contributions we find 

+ K^W-tfl^-^'j/^X!-;) (32) 

where k' — k — (1 — 0)k and 

D lq = k 2 + P(l - (3)Q 2 + m 2 q 



D,„ = (k - k) 2 + (3{l - (3)Q 2 + m 2 



>2q 



k 2 



k' 2 + m 2 

1 + /?(i-/?)Q 2 + g 2 



(33) 



The argument of as is taken as (k +k )+m 2 . We set the quark masses to be m u = = m s = 
and m r — 1.4 GeV. 



3.1 The light quark component of the sea 

We first discuss the calculation of the contribution of the "massless" u, d, s quarks to the 
total sea distribution S. It is necessary to consider three different regions of the k and k' 
integrations of (|3^). 

(a) The contribution from the non-perturbative region k 2 , k' 2 < k 2 is evaluated phenomeno- 
logically assuming that it is dominated by "soft" Pomeron exchange [ 27 1 . The contribution 
is parametrized by the form 

= Sf + Sf + Sf (34) 
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where 

= Sf = 2Sf = Cjpx-° m (l-x) 8 . (35) 

The coefficient Cp is independent of Q 2 (in the large Q 2 region) since the contribution 
arises from the region in which the struck quarks have limited transverse momentum, 

.2 

o- 



(b) In the region k 2 < k 2 , < k' 2 we apply the strong-ordered approximation at the quark-gluon 
vertex and take [28 

Sbox - Si b l(z,k 2 = 0,Q 2 ). (36) 
Then the contribution to (|31|) from this domain becomes 



S ( b) 



f ^s«(^ 2 =o,q 2 )- g r z -x 

J x Z Z \ Z 



(37) 



where the summation over u, d, s is implicitly assumed. The potential collinear singular- 
ities in the on-shell structure function Sb ox are regulated by the cut-off k^. Recall that 

.2 

o- 



K 2 ~ K' 2 > kn. 



(c) In the remaining region, k 2 > k^, eq. (|3^) is left unchanged. To be precise we use the 
perturbative expression for S q (x,Q 2 ). 

3.2 The charm component 

The calculation of the charm component of the sea follows perturbative QCD in all regions. 
To evaluate S q=c we divide the integration over k 2 into the regions k 2 < k 2 . and k 2 > k 2 ,. For 
k 2 < k 2 , which we denote region (b), we use the on-shell approximation to evaluate S^ox- That 
is we calculate Sb ox (z, k 2 = 0, Q 2 ; m 2 ), which is finite due to m c ^ 0. Then (^) gives 

3%{x, Q 2 ) = £ ^ S hox (z, k 2 = 0, Q 2 ; m 2 c ) f ^ / k 2 ) (38) 

where a = (1 + 4m 2 /Q 2 )^ 1 , see (j33|). For k 2 > k 2 ,, which we call region (c), we use the full 
perturbative formula. Thus adding the two contributions we have 



S q = c {x,Q 2 ) = [ — S hox (z, k 2 = 0, Q 2 ] m 2 c )- g (-, k%) 

Jx Z Z \Z J 



where we have used (|l|) which enables S q=c is to be specified in terms of the conventional gluon 
input distribution. 
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3.3 The equation for the quark singlet distribution 

Besides S, the singlet momentum distribution S also contains a valence quark contribution 
V, which is taken from a known set of partons. Thus in summary the singlet distribution is 

S = (S^ + S® + S^) uds + (S® + S®) q=c + V, (40) 

where is phenomenologically parametrized in terms of "soft" Pomeron exchange and the 
terms are determined perturbatively except for the (non-perturbative) input gluon distribution 
at scale k 2 ,. The terms are defined entirely in terms of the unintegrated gluon distribution / 
in the perturbative region k 2 > k 2 ,. Finally, V = x(u va \ + d va \) is the valence quark contribution. 

In order to see the connection with the GLAP evolution of the (light) quark sea we first 
note that 

?Q 2 dS q (x,Q ' 2 ' 
h% dQ'- 

where here S denotes the sum over just the u, d and s quarks. We next recall that the leading 



S q (x,Q 2 ) = S q (x,k 2 ) + j o ^to,2 W 2 , (41) 



twist part of the kx factorization formula (EH), written in the form 



Q2 dS q (x,Q 2 ) = ri dz r dtf Q2 d£%J^, k 2 ,Q 2 ) Jx_ fc2 , ^ 



dQ 2 Jx z J k 2 dQ 2 J \z 

can be reduced to the collinear form [1181 



Q ,a_S^m = asm £ dzPqg(z , asm l g (l, Q ^ (43) 

which incorporates leading hxl/x resummation effects in both the splitting function P qg and in 
the integrated gluon distribution g. Thus ( |4"1~D may be written in the form 

q( n ^ q( f Q2 dQ' 2 a s (Q' 2 ) f l , 
S q (x,Q) = S q {x,k ) + / o — ^— / dz 

(44) 



kl Q' 2 2ir Jx 



P qg (z, a s (Q' 2 )) - z g (Z, Q 12 ) + P„(z, a s (Q' 2 ))S q Q' 2 



where for consistency we have included the S —>■ S contribution to the evolution. This additional 
term is needed to ensure the correct GLAP structure. Of course, at small x we expect S to 
be dominantly driven by the gluon. Equation ([0]) is simply the integral form of the GLAP 
evolution equation for the (light) sea quark (momentum) distribution, S. 

Guided by the GLAP structure, it is clear that we should also add the S — > S contribution 
to the complete equation (flUD based on fcr factorization. Then (|40|) becomes 

£(x,A; 2 ) = S^\x) + J2[ a -Sl x (z,k' 2 = 0,k 2 ] m 2 q )^g(^k 2 > ) +V(x,k 2 ) 

„ Jx Z Z \Z / 



1 



^ dk' 2 a s (k' 2 ) /-i J „ . NO fx ll2 
+ / "775 ^r~~ j dzP qq (z)S uds (- } k' 2 



(45) 



k * fc' 2 2vr Jx ~~~^~'~^\ z - 
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where is given by (|34|) and where the uds subscript indicates that the additional S — > S 
term is only included for the light quarks. This equation for the singlet quark distribution 
E, together with eq. ( [TO] ) for the gluon, form the pair of coupled equations which we solve. 
In this way we can specify the structure function Fi in terms of the parameters of the input 
distributions, and hence determine the values of the parameters by fitting to the data for F2, 
see sections 4-6. 



3.4 k T versus collinear factorization and P qg 

As we have already mentioned, the leading-twist part of the kx factorization formula can 
be rewritten in a collinear factorization form. Once the unintegrated gluon distribution is 
taken as a solution of the BFKL equation and the kx factorization integral is performed over 
the entire available phase-space (i.e. not only over the region corresponding to the strongly 
ordered transverse momenta) then the leading small x effects are automatically resummed in 
the splitting functions and in the coefficient functions. The kx factorization theorem can in 
fact be used as the tool for calculating these quantities p. [| . 

We illustrate this point by using the example of the calculation of the splitting function P qg . 
For simplicity we assume that the coupling as is fixed and that the quarks are massless. We 
begin from the kx factorization formula (f42|) written in moment space 

where we have noted, for massless quarks, that S^ ox is a function of the ratio Q 2 /k 2 . Thus we 
may factorise the convolution over k 2 by taking moments. We find 

^ dS(u;,Q 2 ) 1 /-Woo - - (Q 2 y 



dQ 2 2ni A/2-ioo V k\ 



where S g {u^) and f(u,j) are the Mellin transform of the moment functions S q ox (u , k 2 , Q 2 ) 
and f(uj, k 2 ) i.e. 

S hox (uj, k 2 , Q 2 ) = — d 1 S hox (u, 7) It*) (48) 



2ni Ji/2-ioo \ k 2 

and 



1 /-1/2+ioo k 2 \' 

/>, k 2 ) = — I d7/(w,7) 72 • (49) 



2m Ji/2-ioo \kQ 

Retaining the leading pole 7 99 contribution of f(ui, 7) in the 7 plane, the integrals ( fi"7|) and 
fl49|) can be evaluated to give 



,d S(uj,Q 2 ) 2 ~ ,(Q 2 ^ 99 



Q 2 —^^- = 7g 9 5 box (U, Jgg) C (u , « S ) ( ] (50) 

and 

f(,n k 2 \ = n/..nd,i n*} I 
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/ u2 \ 799 

f(u;,k 2 ) = lgg C(uj,a s ) -j . (51) 



where j gg C is the residue of the pole. The function j gg (as,uj) is the (leading twist) gluon 
anomalous dimension. From (BTj) we see that the integrated gluon is given by 

g(uJ,Q 2 ) = C(uj,a s )r^) (52) 

where g(cj, Q 2 ) is the moment function of the (leading twist part) of the gluon distribution. 
Thus by comparing (|50| ) with the conventional GLAP form 

2 dS q (u, Q 2 ) a s ^ , w „ 2 \ /c-o\ 
Q — dQ 2 = 2^ qg ( ' as)g ^' Q > ( 5 > 

we can identify the moment of the P qg splitting function to be 

Gig — ~ 

— P qg (uj, a s ) = l gg (a s , uj)S^ ok (u, l gg {a Sl u)) (54) 

in the so-called Qq regularization and DIS scheme || which we implicitly adopt. 
In the leading ln(l/a;) approximation we have 

^-P qg (u, as) = (7 BFKL ) 2 ^ ox (" = 0, 7 BFKL ) • (55) 

The anomalous dimension ^/ BFKL has the following expansion [2J| 

which in turn gives for the splitting function P gg 

zP gg (z, a s ) = ^ c n , (57) 

n=l V L L )- 

whereas representation ( |55|) generates the following 

expansion of the splitting function P qg {z, as) at small z 

[asHl/z)] 1 



:P qg (z, as) = zP£\z) + a. £ ^ ^ ^ . (58) 



„=i K n 



The first term on the right hand side vanishes at z = 0. It should be noted that the splitting 
function P qg is formally non- leading at small z when compared with the splitting function P gg . 
For moderately small values of z however, when the first few terms in the expansions (|58l ) and 
(57) dominate, the BFKL effects can be much more important in P qg than in P gg . 

This comes from the fact that all coefficients b n in (|58|) are different from zero, while in 
we have c<i = C3 = ||29|| . The small x resummation effects within the conventional 



QCD evolution formalism have recently been discussed in refs. |3l], 32, 33, 34]. These 
studies already emphasize this point, namely that at the moderately small values of x which 
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are relevant for the HERA measurements, the ln(l/x) resummation effects in the splitting 
function P qg have a much stronger impact on F 2 than do those in the splitting function P gg . 
In particular we should also recall that the BFKL effects in the splitting function P qg can 
significantly affect the extraction of the gluon distribution from the experimental data on the 
slope of the structure function F 2 

Here we also include the subleading ln(l/x) terms which would come from the subleading terms 
in 7 99 etc. Keeping the exact kr factorisation (and not just its large Q 2 limit) we also include 
the non- leading twist contributions to F 2 . They would formally be generated by the contribu- 
tions given by the (non-leading) twist anomalous dimensions. 



4. Numerical analysis and the description of F2 

We now have a closed system of two coupled integral equations for two unknowns. Namely 
equation ([IT]) of Section 2 for the unintegrated gluon distribution f(x, k 2 ) and equation P5| ) of 
Section 3 for the integrated quark singlet (momentum) distribution £(x, k 2 ). The effect of the 
gluon in the perturbative region, k 2 > k 2 , is of special interest. It is the 'dynamo' which drives 
small x physics. 

The advantage of this formulation of the unified BFKL/GLAP equation is that the input 
is well-controlled. We emphasized in Section 2 that the equation for f(x, k 2 ) required only the 
specification of an input form for the integrated gluon, 

xg(x,k 2 ) = N{l-xY, (60) 

say. Moreover, the equation for the singlet S(x, k 2 ) requires as input only the non-perturbative 
sea contribution whose form we assume is given by the "soft" Pomeron 

= Cjpx-° m {l-xf (61) 

and the contributions of fl3~7|) and (|38|) which depend on xg(x, k$) of (|60"D. The choice of 
the exponent -0.08 is motivated by the Regge pomeron intercept found in the analysis of total 



cross section data [[27]]. We choose the exponent of (1 — x) to be 8, typical of the behaviour 
of the sea distribution. In our small x analysis any similar choice would be equally good and 
would not change the quality of the description. 

The valence quark contribution V(x, k 2 ) in (pE5|), which is determined mainly by fixed target 
deep inelastic data, is taken from the leading order GRV set of partons ||. We are there- 
fore able to self-consistently determine f(x,k 2 ) and S(x, k 2 ) as functions of a small number 
of physically motivated parameters. In fact, we have only the two parameters, namely N and 
(3 determining the input gluon distribution flSTjp. The momentum sum rule fixes the value of 



Cjp, which determines the input sea, fl6~T|). The presence of BFKL - like terms means that the 
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momentum sum rule is not exactly conserved. However the violation is small. For example, 
after evolution to Q 2 = 50GeV 2 we find that the sum of the momentum fractions carried by 
the gluon and the light quarks is only increased from 1 to 1.007. We neglect this small violation 
of momentum conservation. 



4.1 The optimum description of the F 2 data at small x 

We determine the values of the input parameters by fitting to the HERA measurements of 
the proton structure function F 2 using 

F 2 = J2e 2 q (S g + V g ), (62) 
which holds in the DIS scheme. We thus have to calculate S q (x,Q 2 ) in terms of the input 



gluon parameters N and /3. To do this we solve the pair of equations ([Tof) and ( 45|) for f(x, k 2 ) 
and k 2 ) using an extension of the method proposed in |35fl . This method incorporates the 
interpolation in two variables x and Q 2 with orthogonal polynomials. Thus the coupled integral 
equations can be transformed into the set of linear algebraic equations and readily solved. In 
this way we can express F 2 (x,Q 2 ) in terms of iV and (3. We then determine the optimum 



)2\ 



values of these parameters by fitting to the HERA [Q] and fixed-target |36j data for F 2 (x,Q 
that are available in the small x domain, x < 0.1. We take a running coupling which satisfies 
a s (M 2 ) = 0.12. 

We actually show the results of two fits. The first is the 'realistic' fit with the kinematic 
constraint imposed (which requires the virtuality of the exchanged gluons along the ladder to 
satisfy \k' 2 \ ~ k'^). Then for comparison we repeat the analysis without imposing the kinematic 
constraint, that is we omit G function in (TO). The quality of the fits are shown in Figs. 2 and 
3, and the parameters given in Table 1. To be precise Figs. 2 and 3 respectively show the 
description of the HI and ZEUS data fl| together with those fixed-target data that occur at 
the same values of Q 2 . 

The fit with the kinematic constraint included (continuous curves) is significantly better 
than that in which it is omitted (shown by the dashed curves). Without the constraint the 
predicted rise of F 2 is a little too steep at the smallest values of x and Q 2 . Over the remainder 
of the x, Q 2 domain the fit (fit 2) gives a good description of F 2 . It is far better, for example, 
than that shown in ref. [|yj . 

The kinematic constraint, which corresponds to subleading ln(l/x) corrections, lowers the 
'hard' pomeron intercept and improves the description of the data, particularly at the smaller 
values of x. In fact the resulting description of F 2 (x,Q 2 ) with just two free parameters (N 
and 0) is excellent, and is comparable, even a little better than, to that achieved in the global 
parton analyses, see, for example, the \ 2 listed in Table 1. Moreover, the overall behaviour of 
the gluon is much more realistic than that of the fit without the kinematic constraint. It gives 



an acceptable description of the WA70 prompt photon data [37], which directly sample the 



gluon at x ~ 0.4. These data were not used to constrain the gluon. For fit 1 the prediction is 
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some 30% above WA70 data, which is within the QCD scale uncertainties, whereas the gluon of 
fit 2 gives a prediction which is a factor of about 2.5 above the data. It is not surprising that the 
gluons are so different in the two fits since they are both contrived to give satisfatory description 
of the measurements of F 2 at small x, despite the fact that the kinematic constraint significantly 
reduces the gluon intercept uo- It is encouraging that it is the description with the kinematic 
constraint that gives the acceptable large x behaviour of the gluon. For completeness we use 
our determination of the unintegrated gluon to compute the conventional gluon distribution 
xg(x, Q 2 ) and compare the result with the gluons of recent sets of partons obtained in GLAP- 
based global analyses of deep inelastic and related data. To be specific the continuous curves 
in Fig. 4 compare the gluon calculated from f(x, k?) of the fit 1 (via eq. (|])) with the gluon 
distributions of the MRS(R2) f|] and GRV || set of partons (shown by the dashed and dotted 
curves respectively). We see that the behaviour of the integrated gluon is very similar to that 
of MRS(R2). This may be expected since the MRS analysis used to the same HERA data 
as those fitted in the present work, whereas these data were not available at the time of the 
GRV analysis. However, we emphasize the different underlying structure of the present analysis 
and the pure GLAP-based descriptions. We shall see below that in the unified BFKL/GLAP 
approach the rise of F 2 is generated essentially by ln(l/x) effects in the off-shell gluon structure, 
F] 9 of Fig. 1. 





Kinematic 


xg = N(l - x)P 




X 2 /datapoint 




constant 


N 13 


Cp 


[393 points] 


Fit 1 


yes 


1.57 2.5 


0.269 


1.07 


Fit 2 


no 


0.85 0.9 


0.269 


1.8 


MRS(R2) 








1.12 



Table 1: The parameters N and j3 determined in the optimum fit to the available data [1], [3(J] 
for F 2 with x < 0.05 and Q 2 > 1.5GeV 2 , without and with the inclusion of the kinematic 
constraint along the gluon ladder. The value of Cjp of ( p5[ ) is also shown, although this is fixed 
in terms of N and (3 by the momentum sum rule. For comparison we also show the x 2 for 
the same set of HERA and fixed-target data obtained in a recent next-to-leading order GLAP 
global parton analysis ||. For both fit 1 and 2 the gluon carries 45% of the proton's momentum 
at the input scale = lGeV 2 . 

4.2 The effect of the ln(l/x) resummation on the gluon 

Fig. 5 shows the behaviour of the unintegrated gluon distribution f(x, k 2 ) as a function of 
k 2 for x = 10 -3 and 10~ 4 . Three different determinations are shown, each of which start from 
the input 

xg(x,k 2 Q ) = 1.57(1 -x) 2 - 5 
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of fit 1 of Table 1. The continuous and dashed curves correspond, respectively, to the behaviour 
with and without the kinematic constraint. The dotted curve is obtained from a GLAP deter- 
mination in which the BFKL kernel in ( |T0D is replaced by the leading order P gg function. That 
is (|IDD is replaced by 



f(x,k 2 



a s (k 2 



2tt 



J x Z \Z 



+ 



dzP gg (z, , 

hi k' 2 



k ' 2 dk' 2 fx , 



r> 



f dzP gq (z)E(-,k' 

J x \Z 



(63) 



where P gg has the usual form 



6 



l-z 



+ z(l-z) + 



l-z) 



+ —5(l-z) 
12 v ; 



N f 

-fS(l-z). 



(64) 



The comparison of the dashed and dotted curves shows that the differences between the BFKL 
and GLAP approaches are not very big, even for the values of x as low as 10 -3 . The differences 
become more obvious when one considers smaller values of x, around 10~ 4 . Even then the 
discrepancies are only visible at lower values of k 2 . This effect can be explained in terms of 
power series expansion in a s /w of the gluon anomalous dimension 



7 bfkl = ^ + C4 ( 



+ ... 



see (|56|) , where 02 = 03 = 05 = 0. The first term of the expansion, which is common to BFKL 
and GLAP, is clearly dominant for the smaller values of as(k 2 ). Thus we confirm the well 
known result that, in the region of moderately small values of x relevant for the HERA mea- 
surements, ln(l/a;) resummation has little effect on the gluon distribution. If the gluon input 
were adjusted to correspond to the optimum fit with the kinematic constraint imposed, then 
the continuous curve would be comparable to the other two. However, a common input is used 
to show the impact of the kinematic constraint. We see that the resulting gluon is smaller and 
less steep. This implies that subleading ln(l/a;) corrections are significant. 



4.3 Effect of ln(l/x) resummation on the structure function F2 

To investigate the various effects of the ln(l/x) terms we compute Fz{x, Q 2 ) using four 
different procedures but with a common input^ 

xg(x,kl) = 1.57(1 -x) 2 - 5 , 

correspoding to fit 1. The four different determinations are shown in Fig. 6 and correspond to 

3 In this calculation we have only included the light quarks u, d, s (which we treat as massless), since we want 
to avoid any dependence on the choice of scale for the heavy quarks. Such a dependence would spoil the clarity 
of the explanation of some effects which we want to stress. Everywhere else we include also the charm quark. 
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(i) The full unified BFKL + GLAP calculation with the kinematic constraint included, 
eqs. (|1(]) and (|45|) , shown continuous curve. 

(ii) Analogous to (i) but without the kinematic constraint (dashed curve). 

(iii) Replace (jlOf) by the pure GLAP equation in the gluon sector, eq. flSlf ), but keep the full 
kx factorization for the quarks (dotted curve). 

(iv) Pure GLAP evolution for both the gluons and the quarks (dot-dashed curve). That is 
instead of (|45|) we use 

S(x, k 2 ) = S {a \x) + S {b) (x, k 2 ) + V(x, k 2 ) + 

where £(fc /2 , k 2 ) is the evolution length and is defined by, 

Z(k' 2 ,k 2 )= r^-asW 2 ). (66) 

Jk' 2 Q 

One can again see from Fig. 6 that the differences between BFKL (with no kinematical con- 
straint) and GLAP evolution in the gluon sector are not very big. The calculations start to 
differ only at x ~ 10 -4 ( dashed and dotted lines). On the other hand when we compare the 
pure GLAP evolution (with the P qg splitting function) with the equations where the entire 
phase space has been taken into account then the differences are much bigger. This implies 
that the leading order terms in aslnl/x present in the gluon off-shell structure function (F^ 9 ) 
are much more important than the terms in the gluon anomalous dimension resulting from the 
BFKL equation. The effect of the kinematic constraint is again evident. It leads to the change 
from the dashed to the continuous curves. Fig. 6 also enables us to see the x values at which 
the effect of the ln(l/a;) summation effects become important. 



4.4 Predictions for F| and Fl 

Once we have determined the parton distributions we can predict the values of other hard 
scattering observables. At small x we see, via the kx factorization theorem, that the observables 
are 'driven' by the unintegrated gluon distribution f(x, k 2 ). Here we calculate F| and F^. 

The charm component F% of F 2 is given by 

F 2 c (x,Q 2 ) = e 2 c S q=c (x,Q 2 ) 
where the charm sea S q=c is calculated from ( p9|) in terms of the gluon. It is the second term on 



the right-hand side of (|39|) which drives the small x behaviour. The predictions are compared 
with the HI measurements [Q of F 2 C in Fig. 7. The percentage of charm in the deep inelastic 
structure function is shown in Fig. 8. At small x we see that F% is an appreciable fraction of 
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F 2 . Recall that in the massless charm limit the fraction would be 0.4, provided that we are 
below the bottom quark threshold. 

The predictions of the longitudinal structure function Fl are shown in Fig. 9. For Fl the 
kr factorization formula can be written in the form 



39 



F L (x,Q' 



71 



\ f * f £) 2 F 2 (y, g») + £ elf * f £) 2 (l - £) M(v , Affl 



1 2 



/(p^J. ( 67 ) 



where the quark box variables D iq and are defined in ([33|). The behaviour of Fi is driven by 
the gluon through the last term. The argument as is taken to be k 2 + k,' 2 + m 2 as before. 



5. Conclusions 

The natural framework for describing 'hard' scattering observables at small x Bjorken x is 
provided by the gluon distribution f(x, kj) unintegrated over its transverse momentum, together 
with the kx factorization theorem. At small x it is only to be expected that the fcr dependence 
should be treated explicitly. 

In the leading ln(l/x) limit, f(x,k^) satisfies the BFKL equation. To make a smooth 
transition to the larger x domain we have studied a modified equation which treats the leading 
ln(l/x) and the leading ln(Q 2 ) terms on an equal footing. Moreover, we arrange the equation 
so that we need only to consider f(x, k?) in the perturbative domain, k\> k 2 . The integrated 
gluon distribution xg(x,ky) is the only non-perturbative input that is required. At small x, 
the singlet quark distribution E(x,Q 2 ) is controlled by f(x,k^) through the g — ^ qq splitting. 
We therefore extend the formalism to a pair of coupled integral equations for / and S which 
embrace both the BFKL leading ln(l/a;) and GLAP leading ln(Q 2 ) contributions in a consistent 
way. A notable feature of the formalism is that we can retain the full perturbative contribution 
of the quark box which contributes to the sea distribution for all Q 2 . In this way we can isolate 
the non-perturbative contribution to a (scaling) sea contribution whose general form is known 
from 'soft' physics. 

An alternative way to unify the BFKL and GLAP formalisms is based on collinear factor- 
ization. It has been shown that we can reduce the (leading twist part of the) BFKL behaviour 
to collinear form in which the the splitting and coefficient functions contain explicit calculable 
series of a™ (asln(l/x)) n terms. This approach has attracted much interest. However, in the 
introduction we stressed the advantage of working with the unintegrated gluon distribution and 
using the k? factorization theorem, and we mentioned some of the limitations of the reduction 
of the BFKL equation to collinear form. Here we simply state some of the points to consider. 
In the 'unintegrated' formalism it is straightforward to identify the perturbative contributions 
which contribute at all Q 2 and so to avoid subsuming them in the input distributions. Second, 
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there is a natural way to introduce the running of as in the BFKL formalism, that has increas- 
ing theoretical support, which for sufficiently small x goes beyond the Renormalisation Group 
behaviour (and so is difficult to implement in the collinear factorization approach). Thirdly, 
the kinematic constraint along the BFKL ladder is easy to implement in the 'unintegrated' 
formalism. Another point is that the BFKL formalism contains all twists, whereas only the 
leading twist is retained in the collinear approach. Last, but not least, the hp factorization ap- 
proach, which we may symbolically write as F 2 = F 2 g <g> /, is easier to implement. The BFKL 
kernel and the off-shell gluon structure function F 2 9 are calculable perturbatively. We simply 
use leading order in as expressions. The ln(l/x) summations are implicit in the integration 
over the entire k\ phase space of the gluon ladder and in the kx factorization integrals. 

We solved numerically the coupled integral equations for f(x, k T ) and E(x, k T ), and we then 
used the kx factorization theorem to calculate F 2 in terms of a two-parameter input form for 
the gluon, xg{x, k^) = N(l — x) 13 . The parameters N and (3 are determined by a fit to the 
available small x data for F 2 . An excellent description is obtained. The data at the smallest 
values of x give support for the presence of the kinematic constraint, as does the extrapolation 
of the gluon to describe the WA70 prompt photon data at x ~ 0.4. Notice that the rise of 
F 2 with decreasing x is purely of perturbative origin in our description, and that we find a 
significant 'BFKL' component. 

The fact that we achieve an excellent two-parameter fit of the small x data for F 2 is not, 
in itself, remarkable. Other equally good phenomenological fits have been obtained. What 
is encouraging is that we have a theoretically well-grounded and consistent formalism which, 
with the minimum of non-perturbative input, is able to give a good perturbative description 
of the observed structure of F 2 . Moreover the BFKL/GLAP components of F 2 are decided 
by dynamics. In this way we have made a determination of the universal gluon distribution, 
f(x,k T ), which can be used, via kx factorization, to predict the behaviour of other small x 
observables. We showed the predictions for F 2 and Fx- ^ 
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Figure captions 



Fig. 1 The diagrammatic representation of the kj- factorization formula = Fj 79 ®/. At lowest 
order in as, the photon-gluon fusion processes (or to be precise the structure functions 
F 79 of the virtual gluon) are given by the quark box shown (together with the crossed 
box). 

Fig. 2 The two-parameter fit to the F2 data at small x using eq. ( pf) for f(x, k 2 ) with (continuous 
curves) and without (dashed curves) the kinematic constraint. The optimum values of 
the parameters N and (3, which describe the input form of the gluon, are given in Table 
1. The figure shows the HI data [JTJ] together with the E665 and NMC measurements [^] 
which occur at the same values of Q 2 . 

Fig. 3 As for Fig. 1, but for the ZEUS measurements [|l[] of F2, together with the E665, NMC 
and BCDMS data [36] which occur at the same values of Q 2 . 



Fig. 4 The continuous curves show the behaviour of the conventional gluon distribution xg(x, Q 2 ) 
corresponding to fit 1, and calculated using eq. (|lj). For comparison we also show the 
gluon distributions of the MRS (R2) [|J (dashed curve) and GRV (dotted curve) sets 
of partons. 

Fig. 5 The unintegrated gluon distribution f(x,k 2 ) as a function of k 2 for x = 10~ 4 and 10~ 3 
obtained by solving the simultaneous equations for f(x,k 2 ) and S(x, k 2 ). The con- 
tinuous and dashed curves are obtained by using the unified BFKL/GLAP equation 
(|T0| ) for f(x,k 2 ) with and without the kinematic constraint respectively. The dotted 



curve corresponds to using GLAP evolution for /, eq. (j63|) . In each case the input 
xg(x,k%) = 1.57(1 — x) 2,5 is used, where k 2 , = lGeV 2 . 

Fig. 6 The light quark contribution to F 2 (x,Q 2 ) for various Q 2 values obtained from solving 
different sets of coupled equations for the gluon / and the quark singlet £ with, in each 
case, the input xg(x, k%) = 1.57(1 — x) 2 ' 5 where k$ = lGeV 2 . The continuous and dashed 
curves come from solving ( |TU| , |4"o] ) with and without the kinematic constraint. The dotted 
curve is obtained using GLAP in the gluon sector, that is (^,^), whereas the dot-dashed 
curve corresponds to pure GLAP evolution, (|63| , |65D . 



Fig. 7 The predictions for F 2 , compared with HI charm data, obtained from the optimum fit 
(fit 1). 

Fig. 8 Ratio F 2 /F 2 for different values of Q 2 obtained from fit 1. 

Fig. 9 The prediction for the structure function Fl as a function of x for different values of Q 2 
using the parameters of fit 1. 
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